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Abstract 

We revisit existence and stability of two-pulse solutions in the fifth-order Korteweg-de Vries 
(KdV) equation with two new results. First, we modify the Petviashvili method of successive it- 
erations for numerical (spectral) approximations of pulses and prove convergence of iterations in a 
neighborhood of two-pulse solutions. Second, we prove structural stability of embedded eigenval- 
ues of negative Krein signature in a linearized KdV equation. Combined with stability analysis in 
Pontryagin spaces, this result completes the proof of spectral stability of the corresponding two- 
pulse solutions. Eigenvalues of the linearized problem are approximated numerically in exponentially 
weighted spaces where embedded eigenvalues are isolated from the continuous spectrum. Approxima- 
tions of eigenvalues and full numerical simulations of the fifth-order KdV equation confirm stability 
of two-pulse solutions related to the minima of the effective interaction potential and instability of 
two-pulse solutions related to the maxima points. 



1 Introduction 

One-pulse solutions (solitons) are commonly met in many nonlinear evolution equations where dispersive 
terms (represented by unbounded differential operators) and nonlinear terms (represented by power 
functions) are taken in a certain balance. Typical examples of such nonlinear evolution equations with 
one-pulse solutions are given by the NLS (nonlinear Schrodinger) equation, the Klein-Gordon (nonlinear 
wave) equation and the KdV (Korteweg-de Vries) equation, as well as their countless generalizations. 

One-pulse solutions are the only stationary (traveling) localized solutions of the simplest nonlin- 
ear evolution equations. However, uniqueness is not a generic property and bound states of spatially 
separated pulses can represent other stationary (traveling) localized solutions of the same evolution 
equation. For instance, two-pulse, three-pulse, and generally iV-pulse solutions exist in nonlinear evolu- 
tion equations with a higher-order dispersion (represented by a higher-order differential operator). The 
prototypical example of such situation is the fifth-order KdV equation in the form, 

ut + Uxxx ~ u xxxxx + 2uu x =0, x GR, t £ R+, (1.1) 

where u:lx M+ i— ► M. and all coefficients of the nonlinear PDE are normalized by a scaling transfor- 
mation. See Bridges & Derks |BD02j for a review of history and applications of the fifth-order KdV 
equation (|l.lj) to magneto- acoustic waves in plasma and capillary- gravity water waves. 

Traveling localized solutions u(x, t) = (f>(x — ct) of the fifth-order KdV equation (|l.lj) satisfies the 
fourth-order ODE 

0( iv ) -0" + c = </> 2 , zeR, (1.2) 
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where z = x — ct is the traveling coordinate and one integration of the fifth-order ODE in z is performed 
subject to zero boundary conditions on 4>{z) and its derivatives as \z\ — > oo. Existence of localized 
solutions (homoclinic orbits) to the fourth-order ODE l|1.2|) was considered by methods of the dynamical 
system theory. See Champneys |C98j for a review of various results on existence of homoclinic orbits in 
the ODE (Ol) . 

In particular, it is proved with the variational method by Buffoni & Sere BS96 and Groves G98 (see 
references to earlier works in C98]) that the fourth-order ODE ()1.2j) has the one-pulse solution cj>(z) for 
c > 0, which is the only localized solution of the ODE (|1.2|) for < c < \ up to the translation <p(z — s) 
for any s G R. The analytical expression for the one-pulse solution is only available for c = < \ 
with 

For c > j, the fourth-order ODE (|1.2j) has infinitely many multi-pulse solutions in addition to the one- 
pulse solution BS96 ( G98 . The multi-pulse solutions look like multiple copies of the one-pulse solutions 
separated by finitely many oscillations close to the zero equilibrium <j) = 0. Stability and evolution of 
multi-pulse solutions are beyond the framework of the fourth-order ODE ljl.2j) and these questions were 
considered by two equivalent theories in the recent past. 

The pioneer work of Gorshkov &; Ostrovsky explains multi-pulse solutions of the fifth-order KdV 
equation from the effective interaction potential computed from the one-pulse solution t G081 ( 

GOP84 . When the interaction potential has an alternating sequence of maxima and minima (which 
corresponds to the case when the one-pulse solution (p(z) has oscillatory decaying tails at infinity), an 
infinite countable sequence of two-pulse solutions emerge with the distance between the pulses near 
the sequence of extremal points. Three-pulse solutions can be constructed as a bi-infinite countable 
sequence of three one-pulse solutions where each pair of two adjacent pulses is located approximately at 
a distance defined by the two-pulse solution. Similarly, iV-pulse solutions can be formed by a (N — 1)- 
infinite countable sequence of N copies of one-pulse solutions. The perturbative procedure in [H081 
has the advantages that both the linear and nonlinear stability of multi-pulse solutions can be predicted 
from analysis of the approximate ODE system derived for distances between the individual pulses. 
Numerical evidences of validity of this procedure in the context of the fifth-order KdV equation are 
reported in |BG97| . 

A different theory was developed by Sandstede |S98j who extended the Lin's work on the Lyapunov- 
Schmidt reductions for nonlinear evolution equations L90 . In this method, a linear superposition of N 
one-pulse solutions 4*(z) = X^=i ~ s j) * s a solution of the ODE (|1.2|) in the case when the distances 
between pulses are infinite (i.e. |sj+i — Sj\ = oo, Vj). The Jacobian of the nonlinear ODE (|1.2j) defines 
a linear self-adjoint operator from i? 2 (R) to L 2 (R): 

H = c-d 2 z +d A z -2(t){z), c>0, (1.4) 

where the unbounded differential part c — d z + d z is positive and bounded away from zero while the 
exponentially decaying potential term — 2<p{z) is a relatively compact perturbation. When 4>(z) is a 
linear superposition of infinitely-separated one-pulse solutions <&(z — Sj), the Jacobian 7i has N zero 
eigenvalues related to the eigenfunctions $'(z — Sj) due to the translational invariance of the ODE H1.2JI . 
The Lyapunov-Schmidt method leads to a system of bifurcation equations for the distances between 
individual pulses. When cj>(z) is the TV-pulse solution with finitely separated pulses (i.e. — Sj\ < oo, 
Vj), one zero eigenvalue of the Jacobian operator TL survives beyond the reductive procedure due to 
the translational invariance of the A^-pulse solution 4>(z), while N — 1 real eigenvalues bifurcate from 
zero. The reduction method gives not only information about existence of multi-pulse solutions but 
also prediction of their spectral stability in the linearized time-evolution problem S98:. The linearized 
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problem for the fifth-order KdV equation takes the form 

d z Hv = Aw, z G M, (1.5) 

where v : M i— ► C is an eigenfunction for a small perturbation of <f>(z) in the reference frame z = x — ct 
and A G C is an eigenvalue. We say that the eigenvalue A is unstable if Re(A) > and v G H 2 (W). 
We say that the eigenvalue A is of negative Krein signature if Re(A) = 0, Im(A) > 0, v G H 2 (U) and 
(Hv,v) < 0. 

Our interest to this well-studied problem is revived by the recent progress in the spectral theory of non- 
self-adjoint operators arising from linearizations of nonlinear evolution equations |CP06| . These opera- 
tors can be defined as self-adjoint operators into Pontryagin space where they have a finite-dimensional 
negative invariant subspace. Two physically relevant problems for the fifth-order KdV equation (jl.lj) 
have been solved recently by using the formalism of operators in Pontryagin spaces. First, convergence 
of the numerical iteration method (known as the Petviashvili's method) for one-pulse solutions of the 
ODE (|1.2j) was proved by the contraction mapping principle in a weighted Hilbert space (which is equiv- 
alent to Pontryagin space with zero index) .I'SO lj. Second, eigenvalues of the spectral stability problem 
in a linearization of the fifth-order KdV equation Ijl.lj) were characterized in Pontryagin space with a 
non-zero index defined by the finite number of neg ative eigenvalues of H |KPn5l IHPOFij . 

Both recent works rise some open problems when the methods are applied to the A^-pulse solutions 
in the fifth-order KdV equation even in the case of two-pulse solutions (N = 2). The successive 

iterations of the Petviashvili's method do not converge for two-pulse solutions. The iterative sequence 
with two pulses leads either to a single pulse or to a spurious solution with two pulses located at an 
arbitrary distance (see Remark 6.5 in |PS04 ). This numerical problem arises due to the presence of 
small and negative eigenvalues of 7i. A modification of the Petviashvili's method is needed to suppress 
these eigenvalues similarly to the work of Demanet & Schlag jQSQS| where the zero eigenvalue associated 
to the translational invariance of the three-dimensional NLS equation is suppressed. We shall present 
the modification of the iterative Petviashvili's method in this article. 

Our numerical method complements the recent work of Yang & Lakoba |YL06| where the steepest 
descent method is developed to suppress any unstable eigenvalues of the iterative method and to enforce 
its convergence. We note that other techniques for numerical approximation of multi-pulse solutions in 
the fifth-order KdV equation have been developed by Champneys and coworkers |CS93| O(.'K03 on the 
basis of a delicate combination of the numerical shooting method and the continuation techniques. 

Open questions also arise when spectral stability of multi-pulse solutions is considered within the 
linear eigenvalue problem (|1.5|) . By either the Gorshkov-Ostrovsky perturbative procedure or the 
Sandstede-Lin reduction method, the small eigenvalues of the Jacobian operator H result in small 
eigenvalues of the linearized operator d z 7i, which are either pairs of real eigenvalues (one of which is 
unstable) or pairs of purely imaginary eigenvalues of negative Krein signature (which are neutrally stable 
but potentially unstable). Both cases are possible in the fifth-order KdV equation in agreement with 
the count of unstable eigenvalues in Pontryagin spaces (see Theorem 6 in .CP06;). (Similar count of 
unstable eigenvalues and eigenvalues of negative Krein signatures was developed for the NLS equations 
in recent papers KKS04, |P05j.) Since the real eigenvalues are isolated from the continuous spectrum 
of the eigenvalue problem (|1.5|) . they are structurally stable and persist with respect to parameter 
continuations. However, the purely imaginary eigenvalues are embedded into the continuous spectrum 
of the eigenvalue problem qi.5J) and their destiny remains unclear within the reduction methods. It is well 
known for the NLS-type and Klein-Gordon- type equations that embedded eigenvalues are structurally 
unstable to the parameter continuations G90 . If a certain Fermi golden rule related to the perturbation 
term is nonzero, the embedded eigenvalues of negative Krein signature bifurcate off the imaginary axis 
to complex eigenvalues inducing instabilities of pulse solutions |CPV05| . (The embedded eigenvalues of 
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positive Krein signature simply disappear upon a generic perturbation CPV05 .) This bifurcation does 
not contradict the count of unstable eigenvalues |KKS04l IP05| and it is indeed observed in numerical 
approximations of various pulse solutions of the coupled NLS equations PY05 . 

From a physical point of view, we would expect that the time evolution of an energetically stable 
superposition of stable one-pulse solutions remains stable. (Stability of one-pulse solutions in the fifth- 
order KdV equation was established with the variational theory |L99j and the multi-symplectic 
Evans function method jBD02t IBDG02 j . ) According to the Gorshkov-Ostrovsky perturbative procedure, 
dynamics of well-separated pulses is represented by the Newton particle's law which describes nonlinear 
stability of oscillations near the minima of the effective interaction potential GOP84 . Therefore, we 
would rather expect (on the contrary to embedded eigenvalues in the linearized NLS and Klein-Gordon 
equations) that the embedded eigenvalues of negative Krein signature are structurally stable in the linear 
eigenvalue problem (|1.5j) and persist beyond the leading order of the perturbative procedure. (Multi- 
pulse solutions of the NLS and Klein-Gordon equations with well-separated individual pulses are always 
linearly stable since the small purely imaginary eigenvalues of the Lyapunov-Schmidt reductions are 
isolated from the continuous spectrum of the corresponding linearized problems YSJOO .) 

Since the count of unstable eigenvalues in CP06 does not allow us to prove structural stability 
of embedded eigenvalues of negative Krein signature, we address this problem separately by using 
different analytical and numerical techniques. In particular, we present an analytical proof of persistence 
(structural stability) of embedded eigenvalues of negative Krein signature in the linearized problem Q1.5JI . 
We also apply the Fourier spectral method and illustrate the linearized stability of the corresponding 
two-pulse solutions numerically. Our analytical and numerical methods are based on the construction 
of exponentially weighted spaces for the linear eigenvalue problem (|1,5[) . (See |PW94j for analysis of 
exponentially weighted spaces in the context of the generalized KdV equation.) Our results complement 
the recent work of Sandstede |S06| and Bridges et al. BCD06 , where two-pulse solutions of the fifth- 
order KdV equation were also addressed. Additionally, Lyapunov-Schmidt reductions of two-pulse 
solutions of the coupled KdV equations were recently considered by Scheel & Wright |HW06 . 

This article is structured as follows. Section 2 contains summary of available results on existence and 
stability of one-pulse and two-pulse solutions of the fifth-order KdV equation (|1.1|) . Section 3 presents 
a modification of the iterative Petviashvili method for convergent numerical approximations of the two- 
pulse solutions in the fourth-order ODE ()1.2j) . Section 4 develops the proof of structural stability of 
embedded eigenvalues in the eigenvalue problem ()1.5j) and numerical approximations of unstable and 
stable eigenvalues in an exponentially weighted space. Section 5 describes full numerical simulations of 
the fifth-order KdV equation (jl.lj) to study nonlinear dynamics of two-pulse solutions. 



2 Review of available results on two-pulse solutions 

Linearization of the ODE (|1.2|) at the critical point (0, 0, 0, 0) leads to the eigenvalues k given by roots 
of the quartic equation, 

k 4 -k 2 + c = 0. (2.1) 

When c < 0, one pair of roots k is purely imaginary and the other pair is purely real. When < c < |, 
two pairs of roots k are real- valued. When c > 4, the four complex- valued roots k are located symmetric 
about the axes. We will use notations ko = Iul(k) > and Ko = Re(ft) > for a complex root of (|2.1|) 
in the first quadrant for c > \. The following two theorems summarize known results on existence of 
one-pulse and two-pulse solutions of the ODE (|l,2j) . 

Theorem 2.1 (One-pulse solutions) Consider the ODE with c > 0. 
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(i) There exists an isolated one-pulse solution cp(z) such that <p £ i? 2 (l)nC 5 (t), <p(— z) = (j)(z), and 
4>(z) — ► exponentially as \z\ — ► oo. Moreover, (j)(z) is C m (M) for any m > 0. 

T/ie Jacobian operator H associated with the one-pulse solution <p{z) has exactly one negative 
eigenvalue with an even eigenfunction and a simple kernel with the odd eigenfunction 4>'(z). 

(Hi) Assume that the map 4>{z) from c > to H 2 (R) is C 1 (R+). Then P'(c) > 0, where P(c) = ||<^||| 2 - 
The linearized operator d z TL has a two-dimensional algebraic kernel in L 2 (M) and no unstable 
eigenvalues with Re(A) > 0. 



Proof, (i) Existence of a symmetric solution (f)(z) in H 2 (M.) follows by the mountain-pass lemma and 
concentration-compactness principle (see Theorem 8 in |G98| and Theorem 2.3 in |L99 ). The equiv- 
alence between weak solutions of the variational theory and strong solutions of the ODE (|1.'2[) is es- 
tablished in Lemma 1 of |G98| and Lemma 2.4 of |L99j . The exponential decay of (j)(z) follows by the 
Stable Manifold Theorem in Appendix A of |BS96j . Global uniqueness of the symmetric solution 4>{z) 
for < c < \ (module the space translation in z) is proved in |AT92j . Local uniqueness of the symmet- 
ric one-pulse solution <j)(z) for c > \ follows from the fact that the Jacobian operator 7i has a simple 
kernel due to the translational invariance (statement (h)). Finally, the smoothness of the function (j)(z) 
is proved from the ODE (|1.2|) by the bootstrapping principle |CCP06j . 

(ii) The Jacobian operator 7i coincides with the Hessian of the energy functional under a fixed value 
of the momentum functional BSS87 . Uniqueness of the negative eigenvalue of 7i follows from the fact 
that the one-pulse solution (j>(z) is a ground state (a global minimizer) of the constrained manifold 
(see Proposition 16 in [£98 ). The ground state for operators with even derivatives and symmetric 
potential functions is always given by the even eigenfunction. The kernel of 7i is simple due to the 
duality principle in Theorem 4.1 of .BS96|- If it is not simple, then global two-dimensional stable and 
unstable manifolds coincide and the time for a homoclinic orbit to go from the local unstable manifold 
to the local stable manifold is uniformly bounded. However, it was shown in |B95j that a sequence 
of homoclinic solutions {u n }neN exists such that the time between local manifolds grows linearly in n. 
Hence, the kernel of Ti is simple with the odd eigenfunction <j>'(z) due to the space translation. 

(iii) Smoothness of the map 4>(z) from c > to H 2 (W) is a standard assumption (see Assumption 
5.1 in |L99j ). The positive value of P'(c) and the absence of eigenvalues of d z TL with Re(A) > follow 
from the stability of one-pulse solutions in Main Theorem of BSS87 , Theorems 4.1 and 5.3 of |L99j 
and Theorem 8.1 of BD02 . Indeed, the only negative eigenvalue of Tt becomes a positive eigenvalue 
in the constrained subspace where the corresponding eigenfunction is orthogonal to (f>(z) |SS90| . The 
two-dimensional algebraic kernel of d z TC follows from the derivatives of the ODE (|1.2jl in z and c: 

H<t>\z) = 0, Hd c <J)(z) = -4>{z). (2.2) 

The algebraic kernel of d z TL is exactly two-dimensional under the condition P'(c) ^ |PW92j . □ 



Theorem 2.2 (Two-pulse solutions) Consider the ODE \l.z}) with c > j. There exists an infinite 
countable set of two-pulse solutions (f>(z) such that (j) £ H 2 {W) n C 5 (]R), (f>(—z) = (f>(z), 4>(z) — > 
exponentially as \x\ — ► oo, and <j)(z) resembles two copies of the one-pulse solutions described in Theorem 
\2.1\ which are separated by small- amplitude oscillatory tails. The members of the set are distinguished 
by the distance L between individual pulses which takes the discrete values {L n } n( z?q. Moreover, for any 
small 5 > there exists 7 > such that 



2im 
ko 



7 



n £ N. 



(2.3) 
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Proof. Existence of an infinite sequence of geometrically distinct two-pulse solutions with the distances 
distributed by (|2.3f) follows by the variational theory in Theorem 1.1 of BS96, under the assumption 
that the single-pulse solution (j>(z) is isolated (up to the space translations). This assumption is satisfied 
by Theorem I2.1f ii). □ 

The following theorem describes an asymptotic construction of the two-pulse solutions, which is used 
in the rest of our paper. 

Theorem 2.3 (Lyapunov Schmidt reductions for two-pulse solutions) Let O \ and let $(z) 
denote the one-pulse solution described by Theorem \2.1\ Let L = 2s be the distance between two copies 
of the one-pulse solutions of the ODE \1.2\) in the decomposition 

<t>(z) = §{z- s) + ${z + s)+ip{z), (2.4) 

where (f(x) is a remainder term. Let W(L) be C m (M + ) function for any m > defined by 

W{L)= [ $ 2 (z)$(z + L)dz. (2.5) 

There exists an infinite countable set of extrema ofW(L), which is denoted by {L n } ne fq. 

(i) Assume that W"{L n ) ^ for a given n 6 N. There exists a unique symmetric two-pulse solution 
cj>(z) = 4> n {z) described by Theorem \2.2\ such that the distance L between individual pulses satisfies 
the bound 

\L-L n \ <C n e- KoL (2.6) 

for some C n > 0. 

(ii) The Jacobian H associated with the n-th two-pulse solution <f> n (z) has exactly two finite negative 
eigenvalues with even and odd eigenf unctions, a simple kernel with the odd eigenfunction 4>' n {z) 
and a small eigenvalue \i with an even eigenfunction, such that 



, 2W"{L n ) 



Q(c) 



< C n e- 2KoL (2.7) 



for some C n > 0, where Q(c) = ||$'||?2 > 0. In particular, the small eigenvalue fi is negative 
when W"(L n ) > and positive when W"(L n ) < 0. 

(Hi) There exists a pair of small eigenvalues A of the linearized operator d z 1i associated with the n-th 
two-pulse solution 4> n {z), such that 



a2 + W"{L n ) 



P'(c) 



< C n e- 2KoL , 



for some C n > 0, where P(c) = H^H^ and P'(c) > 0. In particular, the pair is real when 
W"(L n ) < and purely imaginary (up to the leading order) with negative Krein signature when 
W"(L n )>0. 



Proof. When the tails of the one-pulse solution $(z) are decaying and oscillatory (i.e. when c > ^), the 
smooth function W(L) in (|2.5|) is decaying and oscillatory in L and an infinite set of extrema {L n }neN 
exists. Let us pick L n for a fixed value of n € N such that W'(L n ) = and W"{L n ) ^ 0. 
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(i) When the decomposition (|2.4|) is substituted into the ODE (|1.2j) . one can find the ODE for p(z): 

(c-d 2 + dj- 2<S>(z - s) - 2$(z + s))ip-ip 2 = 2${z - s)<S>{z + s). (2.9) 

Let e = 0{e~ K ° L ) be a small parameter that measures the L°°-norm of the overlapping term <I>(z — 
s)&(z + s) in the sense that for each e > there exist constants Co > and sq > such that 

\\$(z - s)$(z + s)\\ L °o < C e Vs>s . (2.10) 

Denote L = 2s and e*$>(z,L) = 2$(z)$(z + L) and rewrite the ODE ((23|) for (p(z) = ip(z + s): 

(c-d 2 + dj- 2$(z)) <p - 2$(z + L)cp - tp 2 = eV(z, L). (2.11) 

The ODE (|2.11|) is defined in function space H 2 (W), where the Jacobian for the one-pulse solution 

H = c - d 2 + dt - 2*(«) 

has a simple kernel with the odd eigenfunction $\z) by Theorem I2.1f ii). By the Lyapunov-Schmidt 
reduction method (see |GS85| ). there exists a unique solution (p = <p(z, L, e) € H 2 (M.)\ker(Tt), such that 
ip(z, L,0) = and <p(z, L, e) is smooth in e, provided L solves the bifurcation equation F(L, e) = 0, 
where 

F(L,e) = e(&{z),*{z,L)) +2(&{z),${z + L)<p{z)) + {&(z),p 2 (z)) 

- - {<S> 2 (z),d z $(z + L))-e (d L *(z, L),<p{z)) - e (*(z, L),d z (p(z)) + (&(z),p 2 (z)p.l2) 

Since (p(z, L, e) is smooth in e with <p(z, L, 0) = 0, we obtain the expansion 

F(L,e) = -W'(L) + F(L,e), 

where the first term has the order of O(e) and the second term has the order of 0(e 2 ). The state- 
ment follows by the Implicit Function Theorem applied to the scalar equation F(L,e) = under the 
assumption that the root L n of W'(L) is simple. 

(ii) The Jacobian 7i associated with the n-th two-pulse solution <j)(z) = <p n (z) given by the decompo- 
sition (|2.4|) has the form: 

H = c - d 2 + d\ - 2<&(z - s) - 2§{z + s) - 2tp(z) 

In the limit of e = 0, where e = 0(e~ K ° L ), the Jacobian 7i has a double negative eigenvalue and a 
double zero eigenvalue. By a linear combination of eigenfunctions, one can construct one even and 
one odd eigenfunctions for each of the double eigenvalues. By continuity of eigenvalues of self-adjoint 
operators, the double negative eigenvalue splits but the two negative eigenvalues remain negative if e 
is small. By reversibility of the system, eigenfunctions for simple eigenvalues are either even or odd 
and by continuity of eigenfunctions, there is exactly one even and one odd eigenfunctions for the two 
negative eigenvalues. By the translation invariance, the double zero eigenvalue splits into a simple zero 
eigenvalue which corresponds to the odd eigenfunction (ft' n (z) and to a small non-zero eigenvalue of the 
order of O(e) which corresponds to an even eigenfunction. The splitting of the double zero eigenvalue 
in the problem 7iv = [xv is considered by the perturbation theory, 

v(z) = ai*'(2 -s)+ a 2 &{z + s) + V(z), (2.13) 

where (aj, ol^) are coordinates of the projections to the kernel of TL at e = and V(z) is the remainder 
term of the order of O(e). By projecting the eigenvalue problem Tiv = fiv to the kernel of H and 
neglecting the terms of the order of 0(e 2 ), one can obtain a reduced eigenvalue problem: 

VQ{c)ai = -Wax + W"(L n )a 2 , nQ(c)a 2 = W" (L n )a x - Wa 2 , 
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where Q(c) = ||<E>' 




> 0, W"(L n ) is computed from ()2.5j) and 



W = 2 - s)] 2 , + ${z + s))=2 {[^'{z)} 2 , <p{z) + $(z + L)) . 



Since one eigenvalue must be zero with the odd eigenfunction (ft' n (z), the zero eigenvalue corresponds to 
the eigenfunction (|2.13|) with a\ = a<i up to the leading order. By looking at the linear system, we find 
that the zero eigenvalue corresponding to a.\ = 02 exists only if W = W"(L n ). The other eigenvalue at 
the leading order is fi = —2W"{L n )/Q(c). The non-zero eigenvalue of the reduced eigenvalue problem 
corresponds to the even eigenfunction ()2.1M)) with a% = —02- The eigenvalue fi has the order of 0(e). 
By continuity of isolated eigenvalues 7i with respect to perturbation terms of the Lyapunov-Schmidt 
reductions, we obtain the result l|2.7|) . 

(hi) In the limit of e = 0, where e = 0(e~ K ° L ), the linearized operator d z 7i for the n-th two-pulse 
solution 4> n (z) has a four-dimensional algebraic kernel according to the two-dimensional kernel of the 
one-pulse solution (|2.2[) . By the translation invariance, the two-dimensional algebraic kernel survives for 
small e with the eigenfunctions {4>' n (z), d c 4> n (z)}. Two eigenvalues A of the operator d z 7i may bifurcate 
from the zero eigenvalue at the order of 0(y/e). The splitting of the zero eigenvalue in the problem 
dzTLv = Xv is considered by the perturbation theory, 



where («i, 012, Pi, pi) are coordinates of the projections to the algebraic kernel of d z Tt at e = and 
V{z) is the remainder term of the order of 0(e). By projecting the eigenvalue problem d z 7iv = Xv to 
the algebraic kernel of the adjoint operator —7id z and neglecting the terms of the order of 0(Ae, e 2 ), we 
immediately find at the order of 0(A) that Pj = Xctj, j = 1, 2, while (ai, 02) satisfy a reduced eigenvalue 
problem at the order of 0(A 2 , e): 



Isolated eigenvalues d z TC are continuous with respect to perturbation terms of the Lyapunov-Schmidt 
reductions, so that we immediately obtain the result 1)2. 8|) for A £ R when W"(L n ) < 0. In order to 
prove 1)2. 8|l for A G iM when W"(L n ) > 0, we compute asymptotically the energy quadratic form 



where v(z) is given by the eigenfunction 1)2.14)1 with a\ = —02 = 1 and pj = Aay, j = 1, 2. When A G iR 
and W"(L n ) > 0, we have (Ti.v,v) < up to the leading order, such that A is an eigenvalue of negative 
Krein signature. Continuity of the eigenvalues of negative Krein signature (even although the eigenvalues 
A S iM are embedded into the continuous spectrum of d z TC) follows from Pontryagin's Invariant Subspace 
Theorem (Theorem 1 in (UP06 ). By using the exponentially weighted spaces |PW94| . continuity of 
eigenvalues A £ iM. holds and the bound 1)2.8)1 is obtained for an eigenvalue A £ C in a neighborhood of 
the point A £ iM. Note that the continuity of one pair of eigenvalues A £ C satisfying ()2.8)) does not 
exclude the possibility of existence of another pair of eigenvalues A £ C (in the exponentially weighted 
space with a negative weight) which also satisfies 1)2.8)1 . □ 

Remark 2.4 Theorem 12 .Ml is a modification of Theorems 1 and 2 in |S98j (see also L90 ) in applications 
to the ODE ()1.2)) and the spectral problem 1)1.5)1 . The proof is lighten compared to the original proofs 



ai*'(z - s) - a 2 &(z + s) + pid c $(z - s) + p 2 d c <S>(z + s) + V(z) 



(2.14) 




(Hv,v) 



4W"(L n )-P'(c)\X\ 2 + 0(Xe,e 2 ), 
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in |L901 [598 . We note that the persistence of purely imaginary eigenvalues (|2,8|) in the full problem 
(J1.5|) cannot be proved with the Lyapunov-Schmidt reduction method since the essential spectrum of 
d z 7i occurs on the imaginary axis (contrary to the standard assumption of Theorem 2 in |S98| that the 
essential spectrum occurs in the left half-plane.) 

Remark 2.5 The following claim from the Gorshkov-Ostrovsky perturbative procedure |G081llG0P84| 
illustrates the role of W{L) as the effective interaction potential for the slow dynamics of a two-pulse 
solution: 

Claim: Let C\,C2 be some positive constants. For the initial time interval < t < Cie K ° L / 2 and up to 
the leading order 0(e~ K ° L ), the two-pulse solutions of the fifth-order KdV equation 11.1)1 can be written 
as the decomposition 

u(x, t) = $(z -ct- s(t)) + $0 -ct + s(t)) + U(x, t), 

where \\U\\l,cx> < C2e~ K ° L and the slow dynamics of L(t) = 2s(t) is represented by the Newton's particle 
law: 

P'(c)L = -AW'(L). (2.15) 

Although rigorous bounds on the time interval and the truncation error of the Newton's particle law 
were recently found in the context of NLS solitons in external potentials (see }JFG506] and references 
therein), the above claim was not proved yet in the context of two-pulse solutions of the fifth-order 
KdV equation (|l.lj) . We note that perturbation analysis that leads to the Newton's particle law (|2.15|) 
cannot be used to claim persistence and topological equivalence of dynamics of the Newton's particle 
to the full dynamics of two-pulse solutions. 

According to Theorem 12.31 an infinite set of extrema of W(L) generates a sequence of equilibrium 
configurations for the two-pulse solutions in Theorem 12.21 Since P'(c) > by Theorem I2.1f iii). the 
maxima points of W(L) correspond to a pair of real eigenvalues A of the spectral problem Q1.5|) . while 
the minima points of W(L) correspond to a pair of purely imaginary eigenvalues A. The two-pulse 
solutions at the maxima points are thus expected to be linearly and nonlinearly unstable. The two- 
pulse solutions at the minima points are stable within the leading-order approximation (|2.8[) and within 
the Newton's particle law (|2.15j) (a particle with the coordinate L(t) performs a periodic oscillation in 
the potential well). Correspondence of these predictions to the original PDE 1)1.1(1 is a subject of the 
present article. We will compute the interaction potential W(L) and the sequence of its extrema points 
{L n } ng N, as well as the numerical approximations of the two-pulse solutions of the ODE (|1.2|) and of 
the eigenvalues of the operator d z 7i in (|1.5j) . 

3 Iterations of the Petviashvili's method for two-pulse solutions 

We address the Petviashvili's iteration method for numerical approximations of solutions of the fourth- 
order ODE (J1.2|) with c > 0. See review of literature on the Petviashvili's method in PS04]. By using 
the standard Fourier transform in L 1 (M) n L 2 (R) 

0(jfe) = f 0(z)e~ ikz dz, keR, 
Jm. 

we reformulate the ODE (|1.2|) as a fixed-point problem in the Sobolev space H 2 (R): 

m = rjjlvy * e R, (3.1) 
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where 4> 2 (k) can be represented by the convolution integral of (j){k) to itself. An even real- valued solution 
<j)(— z) = 4>{z) of the ODE 1)1.2)1 in H 2 (M.) is equivalent to the even real-valued solution <p{—k) = 4>{k) 
of the fixed-point problem ()3.1|) . Let us denote the space of all even functions in H 2 (JL) by H 2 V (W) and 
consider solutions of the fixed-point problem 1)3.1 Jl in H 2 rir ] 



ev V-" v 

2 



Let }^L be a sequence of Fourier transforms in H 2 V (M) defined recursively by 



where uo(k) € H 2 V (W) is a starting approximation and M n = M[u n ] is the Petviashvili's stabilizing 
factor defined by 

r , Uc + k 2 + k A ) [u{k)fdk , , 

M[u] = ^ y n . (3.3) 

j R u(k)u 2 (k)dk 

It follows from the fixed-point problem ()3.1|) that M[if>] = 1 for any solution <j> £ H 2 V (R). The following 
theorem was proved in (PSOi] and reviewed in [DS06 . 

Theorem 3.1 Let (j){k) be a solution of the fixed-point problem in H 2 V (R). Let Ti. be the Jacobian 
operator evaluated at the corresponding solution 4>{z) of the ODE If H has exactly one 

negative eigenvalue and a simple zero eigenvalue and if 



either 6(z) > or 



inf (b(z) 



< % (3.4) 



then there exists an open neighborhood of <j> in H 2 V (M>), in which <j) is the unique fixed point and the 
sequence of iterations {u n (k)}^ =0 in V3.ty) ~ \3. 5)) converges to <f>. 

Proof. We review the basic steps of the proof, which is based on the contraction mapping principle in 
a local neighborhood of <f> in H 2 v (Nl). The linearization of the iteration map l)3.2|) at the solution <p is 
rewritten in the physical space z£las follows: 

v n+ i(z) = -2a n (j)(z) + v n {z) -{c-d 2 z + d A )~ l Hv n {z) , (3.5) 

where a n is a projection of v n onto (f) 2 in L 2 (M): 

_ {4> 2 ,v n ) 
an ~ (^,0)' 

such that u n = <f> + v n and M n = 1 — a n to the linear order. The operator T = (c — d 2 + d^)~ l, H is a 
self-adjoint operator in Pontryagin space Ho defined by the inner product 

V/,9en : [f, g ] = (( c -d 2 z +dt)f,g). 

See |CP06| for review of Pontryagin spaces and the Pontryagin Invariant Subspace Theorem. Since 
c > 0, the Pontryagin space Ho has zero index and, by the Pontryagin Theorem, the operator T in 
LTo has exactly one negative eigenvalue, a simple kernel and infinitely many positive eigenvalues. The 
eigenfunctions for the negative and zero eigenvalues are known exactly as 

Tcj) = -0, T(j>'(z) = 0. 

Due to orthogonality of the eigenfunctions in the Pontryagin space LTo an d the relation 
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we observe that a n is a projection of v n to in Tig, which satisfies the trivial iteration map: 



"n.+i = 0, n > 1. 



Projection of t> n to 0' in IIo is zero since v n S i^g V (M). As a result, the linearized iteration map (|3.5j) 
defines a contraction map if the maximal positive eigenvalue of T in L 2 (R) is smaller than 2. However, 



<7 T 



L 2 



J - 1 < -2 inf («, (c - d 2 z + 

/ ll«ILa=l 



(3.6) 



If <f>(z) > on z £ M., the right-hand-side of (|3.fi[) is zero. Otherwise, the right-hand-side of IJM.fi|) is 
bounded from above by | |inf ze ]R 0(z)|, which leads to the condition (|3.4|) . □ 



Corollary 3.2 Lei </>(z) 6e a one-pulse solution of the ODE with c > defined by Theorem \2.l\ 
Then, the iteration method IS. 2\ )- ^3~S^) converges to 4>(z) in a local neighborhood of cj) in iT| v (R) provided 
that the condition {3.$ is met. 

The condition (|3.4() is satisfied for the positive exact solution (|1.3|) for c = Since the one-pulse 
solution is positive definite for < c < j |AT92j . it is also satisfied for all values of c £ (0, \). However, 
the solution is sign- indefinite for c > \ , such that the condition Q3.4JI must be checked a posteriori, after 
a numerical approximation of the solution is obtained. 

Besides the convergence criterion described in Theorem 13.11 there are additional factors in the nu- 
merical approximation of the one-pulse solution of the ODE (jl.2|) which comes from the discretization 
of the Fourier transform, truncation of the resulting Fourier series, and termination of iterations within 
the given tolerance bound. These three numerical factors are accounted by three numerical parameters: 



(i) d - the half-period of the computational interval z S [—d, d] where the solution <f>{z) is represented 
by the Fourier series for periodic functions; 

(ii) N - the number of terms in the partial sum for the truncated Fourier series such that the grid size 
h of the discretization is h = 2d/N; 

(iii) e - the small tolerance distance that measures deviation of M n from 1 and the distance between 
two successive approximations, such that the method can be terminated at the iteration n if 

Em = \M n - 1| < £ and Eoo = \\u n+ i - u n ||i°° < e. 

and (j) = u n (z) can be taken as the numerical approximation of the solution 4>(z). 



The numerical approximation depends weakly of the three numerical parameters, provided (i) d is 
much larger than the half-width of the one-pulse solution, (ii) N is sufficiently large for convergence 
of the Fourier series, and (iii) e is sufficiently small above the level of the round-off error. Indeed, the 
constraint (i) ensures that the truncation error is exponentially small when the one-pulse solution is 
replaced by the periodic sequence of one-pulse solutions in the trigonometric approximation SSOOj. The 
constraint (ii) ensures that the remainder of the Fourier partial sum is smaller than any inverse power 
of iV (by Theorem I2.1f i). all derivatives of the function (j)(z) are continuous) T00|. The constraint (iii) 
specifies the level of accuracy achieved when the iterations of the method (|3,2j) - (|3,3[) are terminated. 
While we do not proceed with formal analysis of the three numerical factors (see |DS06j for an example 
of this analysis), we illustrate the weak dependence of three numerical factors on the example of the 
numerical approximation (f>(z) of the exact one-pulse solution ()1.3|) . which exists for c = Mx. Numerical 
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Figure 1: The distance E = \\(j) — <j)\\L°° for the ODE (|l,2j) with c = ^ versus the half-period d of the 
computational interval, the step size h of the discretization, and the tolerance bound e. 

implementation of the iteration method ()3.2|) - ()3.3|) was performed in MATLAB according to a standard 
toolbox of the spectral methods TOO . 

Figure 1 displays the distance E = ||</> — 0||l°° versus the three numerical factors d, h, and e described 
above. The left panel shows that the error E converges to the numerical zero, which is O(10 -15 ) in 
MATLAB under the Windows platform, when the step size h is reduced, while d = 50 and e = 10 -15 
are fixed. The middle panel computed for h = 1 and e = 10 -15 shows that the error E converges to the 
level O(10~ 13 ) when the half-width d is enlarged. The numerical zero is not reached in this case because 
the step size h is not sufficiently small. The right panel computed for h = 1 and d = 50 shows that the 
error E converges to the same level O(10 -13 ) as the tolerance bound e is reduced. In all approximations 
that follow, we will specify h = 0.01, d = 50 and e = 10 -15 to ensure that the error of the iteration 
method (|3.2[) - H3.3j) for one-pulse solutions is on the level of the numerical zero O(10~ 15 ). 

Figure 2 (left) shows the numerical approximation of the one-pulse solutions for c = 4, where the 
small-amplitude oscillations of the exponentially decaying tail are visible. We check a posteriori the 
condition ()3.4|) for non-positive one-pulse solutions |inf^ S ]R 4>(z)\ < 2 for c = 4. Figure 2 (right) displays 
convergence of the errors Em = \M n — 1| and E^ = ||u n +i — ^n||L°° computed dynamically at each 
n as n increases. We can see that the error Em converges to zero much faster than the error E^, 
in agreement with the decomposition of the linearized iterative map (jHHj) into the one-dimensional 
projection a n and the infinite-dimensional orthogonal compliment (see the proof of Theorem 13. II) . In all 
further approximations, we will use the error for termination of iterations and detecting its minimal 
values since is more sensitive compared to Em- 
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Figure 2: One-pulse solutions of the ODE (|1.2|) with c = 4 (left) and convergence of the errors Em and 
to zero versus the number of iterations n. 

Since the numerical approximations <fi(z) of one-pulse solutions can be computed for any value of 
c > 0, one can use <j>(z) for a given c and compute the effective interaction potential (|2.5j) . which defines 
the extremal values {L n } n ^. Theorem 12 . 31 guarantees that the two-pulse solution 4> n (z) consists of two 
copies of the one-pulse solutions separated by the distance L near the point L n where W'(L n ) = and 
W"(L n ) / 0. Table 1 shows the first four values of the sequence {L n }^ =1 for c = 1 (where s n = L n /2 
is the half-distance between the pulses). It also shows the corresponding values from the first four 
numerical approximations of two-pulse solutions <f> n (z) (obtained below) and the computational error 
computed from the difference of the two numerical approximations. We can see that the error decreases 
for larger indices n in the sequence {L n } n£ N since the Lyapunov-Schmidt reductions of Theorem 12.31 
become more and more accurate in this limit. 



two-pulse solution 


effective potential 


root finding 


error 


s = s% 

S = S2 
S = S3 
S = Si 


5.058733328146916 
8.196800619090793 
11.338414567609066 
14.479997655627219 


5.079717398028492 
8.196620796452045 
11.338406246900558 
14.479996635578457 


0.02098406988158 
1.798226387474955 • 10~ 4 
8.320708507980612 • 10~ 6 
1.020048761901649 • 10~ 6 



Table 1: The first four members of the sequence of two-pulse solutions of the ODE 1)1.2(1 with c = 1. 



By Theorem 12 ^I n). the Jacobian operator TC associated with a two-pulse solution <j)(z) has one finite 
negative eigenvalue in the space of even functions H^ V (M.) and one small eigenvalue which is either 
negative or positive depending on the sign of W"(L n ). This small eigenvalue leads to either weak 
divergence or weak convergence of the Petviashvili method in a local neighborhood of (j) m H^ V (M.). 
Even if the small eigenvalue is positive and the algorithm is weakly convergent, the truncation error 
from the numerical discretization may push the small eigenvalue to a negative value and lead thus to 
weak divergence of the iterations. 
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Figure 3 illustrates typical behaviors of the errors Em and E^ versus n for the starting approximation 



u (z) = U (z - s) + U (z + s) 



(3.7) 



where Uq(z) is a starting approximation of a sequence {u n (z)} n( zfq which converges to the one-pulse 
solution $(z) and s is a parameter defined near L n /2 for the n-th two-pulse solution 4> n (z). The left 
panel shows iterations for s near s\ and the right panel shows iterations for s near S2- Since W"{L\) > 
and W"(L2) < 0, the two-pulse solution (fii(z) leads to the weak divergence of the iteration method 
(|3.2|) - (|3.3|) . while the two-pulse solution (j>2(z) leads to the weak convergence of the method. 

At the initial stage of iterations, both errors Em and E^ quickly drops to small values, since the 
starting iterations Uq(z^s) converge to the one-pulse solutions ^(z^fs) while the contribution from the 
overlapping tails of <&(z =p s) is negligible. However, at the later stage of iterations, both errors either 
start to grow (the left panel of Figure 3) or stop to decrease (the right panel). As it is explained above, 
this phenomenon is related to the presence of zero eigenvalue of TL in H^ V (M) which bifurcates to either 
positive or negative values due to overlapping tails of $>(z =F s) and due to the truncation error. At the 
final stage of iterations on the left panel of Figure 3, the numerical approximation u n (z) converges to 
the one-pulse solution <&(z) centered at z = and both errors quickly drop to the numerical zero, which 
occurs similarly to the right panel of Figure 2. No transformation of the solution shape occurs for large 
n on the right panel of Figure 3. 

The following theorem defines an effective numerical algorithm, which enables us to compute the 
two-pulse solutions from the weakly divergent iterations of the Petviashvili's method (|3.2|) - (|3,3|) . 

Theorem 3.3 Let (j> n {z) be the n-th two-pulse solution of the ODE 11.2]) defined by Theorems 1 2. 2\ and 
\2.0A There exists s = s* near s = L n /2 such that the iteration method ^3.2]) - k3.3\) with the starting 
approximation uq{z) = <&{z — s) + <&{z + s) converges to 4> n (z) in a local neighborhood of 4> n in H^ V (R). 

Proof. The iteration operator ()3.2j) - (|3,3|) in the neighborhood of the two-pulse solution <p n in H^ V (R) 
can be represented into an abstract form 



where L{e) has a unit eigenvalue at e = 0, N(v n , e) is C k in v n £ H% v with k > 2 such that N(0, 0) = 
D v N(0, 0) = 0, v n is a perturbation of u n to the fixed point (f> n , and e is a small parameter for two-pulse 
solutions defined in Theorem 12.31 By the Center Manifold Reduction for quasi-linear discrete systems 
(Theorem 1 in J03>), there exists a one-dimensional smooth center manifold in a local neighborhood 
of 6 n in iJg V (M). Let £ be a coordinate of the center manifold such that £ G M, £ — corresponds to 
v = 0, and the dynamics on the center manifold is 



where fj,(e) with /z(0) = 1 is an eigenvalue of the linearized operator T(e) at (j) n in Hg V (R) defined 
in Theorem O and /(£„, e) is C k in £ G R with k > 2 and 7(0,0) = <%/(0,0) = 0. Consider the 
one-parameter starting approximation uq(z) = $(z — s) + &(z + s) in a neighborhood of <p n in H^ V (R), 
where s is close to the value s = s n defined in Theorem 12.31 By the time evolution of the hyperbolic 
component of v n (see Lemma 2 in |J03j ). the sequence v n approaches to the center manifold with the 
coordinate £ n . Iterations of £ ra are sign-definite in a neighborhood of £ = 0. Moreover, there exists 
s\ < s n and S2 > s n , such that the sequences {£n(<si)}neN and {£ n (s2)}neN are of opposite signs. By 
smoothness of v n and £ n from parameter s, there exists a root in between s\ < s* < S2 such that 
£ n (s m ) = for all n G N. □ 



v n+1 = L(e)v n + N(v n , e) 



n G N 



Ca+i = M(e)£„ + /(£„,e) 



n G N 
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Remark 3.4 The proof of Theorem 13.31 does not require that the root s* be unique for the one- 
parameter starting approximation uq(z) = §{z — s) + §{z + s). Our numerical computations starting 
with a more general approximation 1|3.7|) show that the root s* is unique in a neighborhood of s n . 

In order to capture the two-pulse solutions according to Theorem 13.31 we compute the minimum of 
the error for different values of s and find numerically a root s = s* of the function 

/(*) = ™i? (-^oo), 

where no is the first iterations after which the value of -EWty increases (in case of the left panel of Figure 
3) or remains unchanged (in case of the right panel of Figure 3). The numerical root s = s* is found by 
using the secant method: 

Sfc- 2 /(sfc_l) - s fe _ 1 /(s fc _ 2 ) , , 

Sk = 77 s 77 \ • ( 3 -°) 

j(Sfc-l) - /(Sfe-2j 

The Petviashvili's method ()3.2j) - (|3,3|) with the starting approximation ()3.7I) where s is close to the root 
s = s* near the point s = s n converges to the two-pulse solution <j) n {z) within the accuracy of the 
round-off error. 

Figure 4 shows the graph of f(s) near the value s = s\ for c = 1. (The graph of f(s) near s = S2 as 
well as other values of s n look similar to Figure 4.) The left panel shows uniqueness of the root, while 
the right panel shows the linear behavior of f(s) near s = s* which indicates that the root is simple. 
Numerical approximations for the first four values of the sequence {s n } n£ N obtained in this root finding 
algorithm are shown in Table 1. We note that the number of iterations of the secant method (|3.8|) 
decreases with larger values of n, such that = 14 for n = 1, = 12 for n = 2, = 10 for n = 3 
and Nh = 9 for n = 4, while the number of iterations of the Petviashili's method for each computation 
does not exceed 100 iterations. 

Figure 5 shows numerical approximations of the two-pulse solutions for c = 1 and c = 4. We can 
see from the right panel that two-pulse solutions with c = 4 resemble the two copies of the one-pulse 
solutions from the left panel of Figure 2, separated by the small- amplitude oscillatory tails. 

We note that comparison of the numerical approximations of two-pulse solutions <j) n {z) (obtained with 
a shooting method) and the asymptotic approximations from the effective potential ()2.5(l is reported in 
BC97 a . Their numerical results (compare Table 1 in |BC97| with our Table 1) show larger deviations 
between the two approximations (since their ODE has a different normalization compared to our (|1.2|) ) 
but the same tendency that the two approximations become identical in the limit of large separation 
distances s n between the two pulses. 

Finally, the three-pulse and multi-pulse solutions of the fixed-point problem ()1.2j) cannot be approxi- 
mated numerically with the use of the Petviashili's method (|3.2j) . In the space of even functions H^ V (M), 
the three-pulse solutions have two finite negative eigenvalues and one small eigenvalue, while the stabi- 
lizing factor of Theorem 13.11 and the root finding algorithm of Theorem 13.31 can only be useful for one 
finite negative eigenvalue and one zero eigenvalue. The additional finite negative eigenvalue introduces 
a strong divergence of the iterative method (|3.2|) which leads to failure of numerical approximations for 
three-pulse solutions. This numerical problem remains open for further analysis. 

4 Eigenvalues of the stability problem for two-pulse solutions 

We address spectral stability of the two-pulse solution <fi n (z) by analyzing the linearized problem (|1.5j) . 
where the operator TL : H 2 (W) i— > L 2 (M) is the Jacobian operator (jl.4[) at cj} n (z). 
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Number of iterations Number of iterations 

Figure 3: Errors Em and E^ versus the number of iterations n for the starting approximation (|3,7[) 
with s = 5.079 (left panel) and s = 8.190 (right panel). The other parameters are: c = 1, d = 50, 
h = 0.01 and s = 10" 15 . 
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Figure 5: Numerical approximation of the first four two-pulse solutions of the ODE (|1,2|) for c = 1 (left) 
and c = 4 (right). 

By Theorem I2.3f ii). operator Ti has two finite negative eigenvalue, a simple kernel and one small 
eigenvalue, which is negative when W"(L n ) > and positive when W"(L n ) < 0. Persistence (struc- 
tural stability) of these isolated eigenvalues beyond the leading order (|2.7|) is a standard property of 
perturbation theory of self-adjoint operators in Hilbert spaces (see Section IV. 3. 5 in |K95j ). 

By Theorem 12 . 'S( m) . operator d z 7i has a pair of small eigenvalues, which are purely imaginary when 
W"(L n ) > and real when W"(L n ) < 0. We first prove that no other eigenvalues may induce instability 
of two-pulse solutions (i.e. no other bifurcations of eigenvalues of d z TC with Re(A) > may occur in 
if 2 (R)). We then prove persistence (structural stability) of the purely imaginary eigenvalues beyond 
the leading order l|2.8|) . Combined together, these two results lead to the theorem on spectral stability 
of the two-pulse solution (fi n (z) which corresponds to L n with W"(L n ) > 0. 

Theorem 4.1 Let -/V rea i be the number of real positive eigenvalues of the linearized problem U.5\) in 
H 2 (M) , N comp be the number of complex eigenvalues in the first open quadrant, and N^ ag be the number 
of simple positive imaginary eigenvalues with (TCv, v) < 0, where v(x) is the corresponding eigenfunction 
for A G Assume that no multiple imaginary eigenvalues exist, the kernel of TC is simple and 

P'(c) > ; where P = \\</)\\ 2 L 2. Then, 

N ie a + 2N oamp + 2Nr as = n{H)-l, (4.1) 

where n(TL) is the number of negative eigenvalues ofTL in H 2 (W). 

Proof. The statement is equivalent to Theorem 6 in _CPQ6j in the case {J-L~ l 4>, (ft) = —(d c cp,(p) = 
— iP'(c) < 0. The result follows from the Pontryagin Theorem in Pontryagin space Tl K , where k = n{TL). 
□ 
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Corollary 4.2 Let <f>(z) = $(z) be a one-pulse solution defined by Theorem \2.1\ Then, it is a spectrally 
stable ground state in the sense that iV rea i = N comp = iV-~ ag = 0. 

Remark 4.3 It is shown in Lemma 4.12 and Remark 4.14 in tCPQ6] that multiple imaginary eigenvalues 
may only occur if (TCv, v) = such that n(TL) > 2 is a necessary condition for existence of multiple 
eigenvalues (with -P'(c) > 0). No multiple imaginary eigenvalues exists for the one-pulse solution &(z). 

Corollary 4.4 Let 4>(z) = (p n (z) be a two-pulse solution defined by Theorem \2.<A Then, 

(i) the solution corresponding to L n with W"(L n ) < is spectrally unstable in the sense that N iea \ = 1 
and -/Vcomp = ^i m ag = f or sufficiently small 0(e~ K ° Ln ) 

(ii) the solution corresponding to L n with W"(L n ) > satisfies N Iegi \ = and N comp + A^~ ag = 1 for 
sufficiently small 0(e~ K ° Ln ) . 

Proof. It follows from Theorems 12. II and 12 .31 for sufficiently small 0(e~ K ° Ln ) that P'(c) > 0, the kernel 
of TL is simple for W"{L n ) ^ 0, and the only pair of imaginary eigenvalues with (TCv, v) < in the case 
W"(L n ) > is simple. Therefore, assumptions of Theorem 14. II are satisfied for the two-pulse solutions 
4> n {z) with W"(L n ) / 0. By the count of Theorem E^ii) , n{TL) = 3 for W"(L n ) > and n(H) = 2 for 
W"(L n ) < 0. Furthermore, persistence (structural stability) of simple real eigenvalues of the operator 
d z TL in H 2 (R) follows from the perturbation theory of isolated eigenvalues of non-self-adjoint operators 
(see Section VIII.2.3 in [EI]). □ 

There exists one uncertainty in Corollary I4.4f ii) since it is not clear if the eigenvalue of negative 
Krein signature in Theorem I2.3f iii) remains imaginary in N-^ or bifurcates to a complex eigenvalue in 
JV comp . This question is important for spectral stability of the corresponding two-pulse solutions since 
the former case implies stability while the latter case implies instability of solutions. We will remove 
the uncertainty and prove that A?;~ ag = 1 and iV C om P = for sufficiently small 0(e _K ° L ). To do so, we 
rewrite the linearized problem (J1.5|) in the exponentially weighted space |PW94j : 



which acts on the eigenfunction v a (z) = e az v(z) G H 2 (JH). The absolute continuous part of the spectrum 
of C a is located at A = X a (k), where 




(4.2) 



The linearized operator d z TL transforms to the form 



C a = (d z - a) (c - (d z - a) 2 + (d z - a) 4 - 2<j>{z)) 



(4.3) 



A Q (fc) = (ik — a)(c — (ik — a) 2 + (ik — a) 4 ) 



feel. 



(4.4) 



A simple analysis shows that 



— Re(A a (/c)) 



2ak(10k 2 - 10a 2 +3), 



c - 3a 2 + 5a 4 + 3k 2 (l - 10a 2 ) + 5k 4 . 



The following lemma gives a precise location of the dispersion relation A = X a (k) on A G C. 
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Lemma 4.5 The dispersion relation A = X a (k) is a simply- connected curve located in the left half-plane 
of A £ C if 

< a < —=, c > -. (4.5) 
x/10 4 



Proof. The mapping A; Im(A a ) is one-to-one provided that c — 3a 2 + 5a 4 > and 1 — 10a 2 > 0. 
Since c — 3a 2 + 5a 4 reaches the minimum value on a £ 0, —t= at the right end a = —k= and the 



10 

minimum value is positive if c > i, the first inequality is satisfied under Q4.5JI . The second inequality 
is obviously satisfied if \a\ < -7=. The mapping k 1— ► Re(A Q ) has a single extremal point at /c = 



provided 3 — 10a 2 > 0, which is satisfied if |a| < -7=. The extremal point is the point of maximum and 



the entire curve is located in the left half-plane ofAeCif0<a< -4= . □ 

v 10 



The following two lemmas postulate properties of eigenfunctions corresponding to embedded eigen- 
values of negative Krein signature. 

Lemma 4.6 Let vq(z) be an eigenfunction of d z TL for a simple eigenvalue Ao £ zR+ in H 2 (Nl). Then, 
Aq 6 is also an eigenvalue in H 2 (R) for sufficiently small a. 



Proof. Let k = ko £ M be the unique real root of the dispersion relation Ao(A;) = Ao (with a = 0) 
for a given eigenvalue Ao £ The other four roots k = ^2,3,4 for a given Ao £ iM+ are complex 

with |Re(/cj)| > kq > 0. By the Stable and Unstable Manifolds Theorem in linearized ODEs |CL55j . 
the decaying eigenfunction vq(z) £ H 2 (JH) is exponentially decaying with the decay rate greater than 
kq > and it does not include the bounded term e tk ° z as z — > ±00. By construction, v a (z) = e az vo(z) 
is also exponentially decaying as z — > ±00 for sufficiently small |a| < kq. Since vo £ H 2 (W) and due to 
the exponential decay of v a (z) as \z\ — > 00, we have u Q £ H 2 (R) for any small a. □ 



Lemma 4.7 Lei vq(z) £ Lf 2 (M) 6e an eigenfunction of d z TL for a simple eigenvalue Ao £ iM + with 
(TIvq^vq) < 0. Then, there exists wo £ Lf 2 (R), such that vq = w' (x) and wo(z) is an eigenfunction of 
7id z for the same eigenvalue Aq. Moreover, (wq,vq) £ iR + . 



Proof. Since TL : H 2 (^l) L 2 (IR), the eigenfunction vq(z) of the eigenvalue problem 8 z TLvq = XqVq for 
any Ao / must satisfy the constraint J R vo(z)dz = 0. Let vo = w' (z). Since vq{z) decays exponentially 
as \z\ —> 00 and (1, Wo) = 0> then wo(z) decays exponentially as \z\ — > 00, so that w £ H 2 (R). By 
construction, TL8 z wq = TLvq = Aq / vo(z)dz = Xwq. The values of (wq, vq) are purely imaginary as 



(w ,v ) = / WQVodz = I w d z w dz = - / w d z w dz = - / w v dz = -(w ,v ). 
Jr Jr jr Jr 

Since TLvq = XqWq with Ao £ iM. + and (Hvo,vq) < 0, we have (wo,vq) = Xq 1 (Hvq,vq) £ iR + . □ 

The following theorem states that the embedded eigenvalues of negative Krein signature are struc- 
turally stable in the linearized problem ()1.5|) . 



Theorem 4.8 Let Ao £ M + be a simple eigenvalue of d z TC with the eigenfunction vo £ i? 2 (lR) such 
that (TCvq,vq) < 0. Then, it is structurally stable to parameter continuations, e.g. for any V £ L°°(R) 
and sufficiently small 5, there exists an eigenvalue X$ £ iM+ of d z (TL + 5V{z)) in H 2 {M), such that 
|A«5 — Aq I < C5 for some C > 0. 
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Proof. By Lemma 14,61 Ao is also an eigenvalue of C a in H 2 (R) for sufficiently small a. Let a be 
fixed in the bound (|4.5I) . There exists a small neighborhood of Ao, which is isolated from the absolute 
continuous part of the spectrum of C a . By the perturbation theory of isolated eigenvalues of non-self- 
adjoint operators (see Section VIII. 2. 3 in |K95j ). there exists a simple eigenvalue A,5 of d z (7i + 5V(z)) 
in H 2 (R) for the same value of a and sufficiently small 5 in a local neighborhood of Ao, such that 
|A<5 - A 1 < C6 for some C > 0. 

It remains to show that the simple eigenvalue Xs is purely imaginary for the same value of a > 0. 
Denote the eigenfunction of d z (7i + 5V(z)) in H 2 (R) for the eigenvalue A^ by v$(z), such that e az vs G 
H 2 (R). If vs $l H 2 (R), then the count of eigenvalues (|4.1j) in H 2 (R) is discontinuous at 5 = 0: the 
eigenvalue Ao in the number at 5 = disappears from the count for 5^0. If vg G H 2 (R), then 

(1,^) = and since v$(z) is exponentially decaying as \z\ — ► oo, there exists w$(z) G H 2 (R) such 
that vs = w's( z )- The 2-form (ws,vs) is invariant with respect to the weight a since if e az vs(z) is an 
eigenfunction of d z (7i + 5V(z)) in // 2 (M) for the eigenvalue Xs (i.e. G H 2 (R)), then e~ a2 w«5(z) is an 
eigenfunction of (7i + 5V(z))d z in H 2 (R) for the same eigenvalue A,5 (i.e. ws G H^^R)). Computing 
(w$,vs) at a = 0, we have 

X s (w s ,v s ) = (Ttvs, vs) el, 

which indicates that (ws,vs) = if Re(A5) ^ and Im^) > 0. Since (wo,vo) ^ by Lemma l4~71 and 
(ws,vs) is continuous in 5, the case Re(A£) ^ is impossible. Therefore, Re(A,5) = is preserved for 
small non-zero 5. □ 

Corollary 4.9 Let (p(z) = (f> n (z) be a two-pulse solution defined by Theorem \2.^ that corresponds to L n 
with W"(L n ) > 0. Then, it is spectrally stable in the sense that N rca \ = N comp = and A^ ag = 1 for 
sufficiently small 0(e~ K ° Ln ). 

Remark 4.10 Using perturbation theory in H 2 (W) for a fixed value a > 0, one cannot a priori exclude 
the shift of eigenvalue Ao to A,5 with Re(A,5) > 0. Even if vq(z) for Ao contains no term e lk ° z as z — ► — oo 
(see Lemma l4.6jl . the eigenfunction 1^(2:) for Xs may contain the term e iksZ as z — > —00 with Im(ks) < 
and lim,5^o ks = ko £ K. However, when Theorem 14.81 holds (that is under the assumptions that 
vq £ ff 2 (M) and (Ttvo,vo) < 0), the eigenvalue Xs remains on iR and the eigenfunction vg(z) must have 
no term e lksZ with ks G M as z — > — cxd for any sufficiently small <5. The hypothetical bifurcation above 
can however occur if dq ^ iif 2 (]R) but dq G iif 2 (lR) with q > 0. We do not know any example of such 
bifurcation in the case vq $ H 2 (R) but (Hv ,v ) > 0. 

Remark 4.11 When the potential is symmetric (i.e. <ft(—z) = 4>(z)), there exists a remarkable symme- 
try of eigenvalues and eigenfunctions of the stability problem d z TCv = Xv: if v(z) is an eigenfunction for 
A, then v{—z) is the eigenfunction for —A. If Ao G iR is a simple eigenvalue and vq G H 2 (R) with a > 0, 
the above symmetry shows that vq G H^ a (R) with —a < 0. If Re(Aa) > and vg G H 2 (R), then —vg(—z) 
is an eigenfunction of the same operator in H^^R) with —a < for eigenvalue Re(— A^) = — Re(A<5) 
and Im(— Xs) = Im(A5). Thus, the hypothetical bifurcation in Eemark 14. 101 implies that the embedded 
eigenvalue Ao G iR may split into two isolated eigenvalues A,5 and — A,5 as 5 7^ 0. Theorem 14.81 shows 
that such splitting is impossible in H 2 (R) if (TLvq, Vq) < 0. 

We confirm results of Corollaries 14.41 and 14.91 with numerical computations of eigenvalues in the 
linearized problem (|1.5|) , Throughout computations, we use the values a = 0.04 and c = 1, which 
satisfy the constraint (|4.5|) . The spectra of the operators TL in H 2 (R) and d z H in H 2 (R) are computed 
by using the Fourier spectral method. This method is an obvious choice since the solution <p(z) is 
obtained by using the spectral approximations in the iterative scheme (|3.2[) ~ ([3.3j) . As in the previous 
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section, we use numerical parameters d = 100, h = 0.01 and e = 10 for the Petviashvili's method 

Eigenvalues of the discretized versions of the operators 7i and C a are obtained with the MATLAB 
eigenvalue solver eig. The spectra are shown on Figure 6 for the two-pulse solution <j)\(z) and on 
Figure 7 for the two-pulse solution 4>2(z). The inserts show zoomed eigenvalues around the origin and 
the dotted line connects eigenvalues of the discretized operators that belong to the absolutely continuous 
part of the spectra. Figures 6 and 7 clearly illustrate that the small eigenvalue of H is negative for 
4>i (z) and positive for 4>2{z)-, while the pair of small eigenvalues of C a is purely imaginary for 4>i{z) and 
purely real for (j)2(z). This result is in agreement with Corollaries 14.41 and 14.91 We have observed the 
same alternation of small eigenvalues for two-pulse solutions 4>3(z) and (p4(z), as well as for other values 
of parameters c and a. 

The numerical discretization based on the Fourier spectral method shifts eigenvalues of the operators 
7i and C a . In order to measure the numerical error introduced by the discretization, we compute the 
numerical value for the " zero" eigenvalue corresponding to the simple kernel of 7i and the double zero 
eigenvalue of C a . Table II shows numerical values for the "zero" and small eigenvalues for two-pulse 
solutions (ftn(z) with n = 1,2,3,4. It is obvious from the numerical data that the small eigenvalues are 
still distinguished (several orders higher) than the numerical approximations for zero eigenvalues for 
n = 1,2,3 but they become comparable for higher-order two-pulse solutions n > 4. This behavior is 
understood from Theorem 12 . 31 since the small eigenvalues becomes exponentially small for larger values 
of s (larger n) in the two-pulse solution (|2.4j) and the exponentially small contribution is negligible 
compared to the numerical error of discretization. 







<h{z) 




4> 4 (z) 


" Zero" eigenvalue of TL 


1.2156 • 10~ y 


2.6678 • 10~ y 


1.4742 • 10~ 9 


1.8941 • 10~ 9 


Small eigenvalue of TL 


1.7845 • 10~ 2 


7.6638 • 10~ 5 


3.3342 • 10"' 


2.9212 • 10~ 9 


" Zero" eigenvalues of C a 


0.3653 • 10" 5 


0.5321 • 10~ 5 


0.7832 • 10~ 5 


1.2374 • 10~ 5 


Re of small eigenvalues of C a 


4.5291 • 10~ e 


3.2845 • 10~ 3 


6.3256 • 10~ 5 


1.6523 • 10~ 5 


Im of small eigenvalues of C a 


0.5021 • 1Q- 1 


1.1523 • 10~ 8 


2.1672 • 10" 4 


5.4443 • 10 _e 



Table II: Numerical approximations of the zero and small eigenvalues of operators TL and C a for 
the first four two-pulse solutions with c = 1, a = 0.04, d = 100, h = 0.01 and e = 10~ 15 . The absolute 
values of real and imaginary eigenvalues are shown. 

We have confirmed numerically the analytical predictions that all two-pulse solutions corresponding 
to the points L n with W"(L n ) < (which are maxima of the effective interaction potential) are unstable 
with a simple real positive eigenvalue, while all two-pulse solutions corresponding to the points L n with 
W"(L n ) > (which are minima of the effective interaction potential) are spectrally stable. The stable 
two-pulse solutions are not however ground states since the corresponding linearized problem has a pair 
of purely imaginary eigenvalues of negative Krein signature. 

5 Nonlinear dynamics of two-pulse solution 

The Newton's particle law H2.15j) is a useful qualitative tool to understand the main results of our article. 
Existence of an infinite countable sequence of two-pulse solution {4> n (z)} ne z is related to existence of 
extremal points {L n } n£ z of the effective potential function W(L), while alternation of stability and 
instability of the two-pulse solutions is related to the alternation of minima and maxima points of 
W(L). It is natural to ask if the Newton's law (|2.15|) extends beyond the existence and spectral stability 
analysis of two-pulse solutions in the fifth-order KdV equation (|1.1|) . In particular, one can ask if the 
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Figure 6: Numerical approximations of the spectra of operators 7i and C a for the two-pulse solution 
4>i(z) with c = 1 and a = 0.04. The insert shows zoom of small eigenvalues and the dotted curve 
connects eigenvalues of the continuous spectrum of C a . 
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purely imaginary (embedded) eigenvalues of the linearized problem (|1.5|) lead to nonlinear asymptotic 
stability of two-pulse solutions or at least to their nonlinear stability in the sense of Lyapunov. From 
a more technical point of view, one can ask whether the Newton's law (|2.15|) serves as the center 
manifold reduction for slow nonlinear dynamics of two-pulse solutions in the PDE and whether 

solutions of the full problem are topologically equivalent to solutions of the Newton's law. While we do 
not attempt to develop mathematical analysis of these questions, we illustrate nonlinear dynamics of 
two-pulse solutions with explicit numerical simulations. 

The numerical pseudo-spectral method for solutions of the fifth-order KdV equation I)1.1J) is described 
in details in M V99j . The main idea of this method is to compute analytically the linear part of the 
PDE by sing the Fourier transform and to compute numerically its nonlinear part by using an 

ODE solver. Let u(k, t) denote the Fourier transform of u(x, t) and rewrite the PDE f| 1 . If) in the Fourier 
domain: 

u t = i(k 3 + k 5 )u - iku 2 . (5.1) 

In order to compute u 2 (k,t) we evaluate u 2 {x,t) on x G M and apply the discrete Fourier transform. 
Substitution u = s(k,t)e l ( k transforms the evolution equation (J5.1|) to the form: 

a t = -ike-^ k3+k5 ^(k, t). (5.2) 

The fourth-order Runge-Kutta method is used to integrate the evolution equation (|5.2|) in time with 
time step At. To avoid large variations of the exponent for large values of k and t, the substitution 
above is updated after m time steps as follows: 

u = s m (k,t)e i{k3+k5){t - mAt \ mAt < t < (m + I) At. (5.3) 

The greatest advantage of this numerical method is that no stability restriction arising from the linear 
part of (|5.1|) is posed on the timestep of numerical integration. On contrast, the standard explicit 
method for the fifth-order KdV equation has a serious limitation on the timestep of the numerical 
integration since the fifth-order derivative term brings stiffness to the evolution problem. The small 
timestep would be an obstacle for the long time integration of the evolution problem due to accumulation 
of computational errors. 

Numerical simulations of the PDE (|5.1[) are started with the initial condition: 

u(x,0) = $(x - s) + 3>(x + s), (5.4) 

where <3?(x) is the one-pulse solution and 2s is the initial separation between the two pulses. The one- 
pulse solution <&(x) is constructed with the iteration method H3.2j) - (|3.3|) for c = 4. The numerical factors 
of the spectral approximation are L = 100, N = 2 12 , e = 10 -15 , while the timestep is set to At = 10 . 

Figure 8 shows six individual simulations of the initial-value problem (|5.1|) and (|5.4|) with s = 2.3, 
s = 2.8, s = 3.6, s = 4.2, s = 4.5 and s = 4.7. Figure 9 brings these six individual simulations on the 
effective phase plane (L, L) computed from the distance L(t) between two local maxima (humps) of the 
two-pulse solutions. 

When the initial distance (s = 2.3) is taken far to the left from the stable equilibrium point (which 
corresponds to the two-pulse solution <pi(x)), the two pulses repel and diverge from each other (trajectory 
1). When the initial distance (s = 2.8) is taken close to the left from the stable equilibrium point, we 
observe small-amplitude oscillations of two pulses relative to each other (trajectory 2). When the 
initial distances (s = 3.6) and (s = 4.2) are taken to the right from the stable equilibrium point, we 
continue observing stable oscillations of larger amplitudes and larger period (trajectories 3 and 4). The 
oscillations are destroyed when the initial distances are taken close to the unstable equilibrium point 
(which corresponds to the two-pulse solution (f>2(x)) from either left (s = 4.5) or right (s = 4.7). In 
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either case, the two pulses repel and diverge from each other (trajectories 5 and 6). Ripples in the 
pictures are due to radiation effect and the numerical integration does not make sense after t « 500 
because the ripples reach the left end of the computational interval and appear from the right end due 
to periodic boundary conditions. 

The numerical simulations of the full PDE problem (|1.1|) indicate the validity of the Newton's particle 
law (|2.15|) . Due to the energy conservation, all equilibrium points in the Newton's law are either 
centers or saddle points and the center points are surrounded by closed periodic orbits in the interior 
of homoclinic loops from the stable and unstable manifolds of the saddle points. Trajectories 2,3, 
and 4 are taken inside the homoclinic orbit from the saddle point corresponding to 4>2(x) and these 
trajectories represent periodic oscillations of two-pulse solutions near the center point corresponding 
to 4>i(x). Trajectories 1 and 6 are taken outside the homoclinic orbit and correspond to unbounded 
dynamics of two-pulse solutions. The only exception from the Newton's law (|2.15JI is trajectory 5, 
which is supposed to occur inside the homoclinic loop but turns out to occur outside the loop. This 
discrepancy can be explained by the fact that the Newton's law (|2.15|) does not exactly represent the 
dynamics of the PDE (|5.1j) generated by the initial condition (|5.4|) but it corresponds to an asymptotic 
solution after the full solution is projected into the discrete and continuous parts and the projection 
equations are truncated (see details in |JFGS06| in the context of the NLS equations). 

Summarizing, we have studied existence, spectral stability and nonlinear dynamics of two-pulse solu- 
tions of the fifth-order KdV equation. We have proved that the two-pulse solutions can be numerically 
approximated by the Petviashili's iterative method supplemented with a root finding algorithm. We 
have also proved structural stability of embedded eigenvalues with negative Krein signature and this 
result completes the proof of spectral stability of two-pulse solutions related to the minima points of 
the effective interaction potential. The validity of the Newton's particle law is illustrated by the full 
numerical simulations of the fifth-order KdV equation which show agreement of slow nonlinear 

dynamics of two-pulse solutions with predictions of the Newton's particle law. 
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right). 
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